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Abstract 

r~ | ■ The three-loop QCD contributions to the vacuum polarization functions of the 

Qh! Z and W bosons at zero momentum are calculated. The top quark is considered 

O ,| to be massive and the other quarks massless. Using these results, we calculate 

the correction to the electroweak p parameter. 

All computations are done in the framework of dimensional regularization 
as well as regularization by dimensional reduction. We use recurrence relations 
^ obtained by the method of integration by parts to reduce all integrals to a small 

set of master integrals. 

A comparison of the two-loop and three-loop QCD corrections to the p pa- 
rameter is performed. 
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Owing to the apparent discovery of the top quark |]J with a mass of 174 ± 10 GeV, 
the prospects grow to test the Standard Model on an even higher level of precision than 
it was possible by now. In particular, the precise knowledge of top-mass effects will 
allow us to obtain better limits on virtual Higgs effects (and thus, indirectly, on the 
Higgs mass), and possibly, on new physics. For this reason, a great deal of work has 
been devoted to the study of the top-mass effects in higher-loop radiative corrections 
of various electroweak parameters. Mostly, in these studies the top mass is assumed 
to be large compared to all other masses so that the latter can be put equal to zero 
from the very beginning (see 0). In || |4[] also the Higgs mass m# was kept as an 
independent parameter, and both limits, niu <C fn t and niu 3> ni t , were studied. 

In the Standard Model there are two different sources of corrections which become 
large (~ G^rri 2 ) in the limit of a heavy top, owing to the large top-bottom mass ratio: 
the Z and W self-energies (affecting, in particular, the p parameter) |J and the Zbb 
vertex 0. Experimentally, these effects are best accessible in e + e _ — > ff near the 
Z resonance measured at LEP1 and the on-resonance asymmetries measured at LEP1 
and the SLAC e + e~ linear collider SLC. 

In the present paper we are concerned with the heavy-top QCD-corrections to the 
electroweak p parameter in three-loop approximation. The p parameter is defined 
as the ratio of the neutral-current to charged-current amplitudes at zero momentum 
transfer: 

Gnc(0) 1 m 
9 G cc (0) 1-Ap' [> 

where the leading fermion contribution to Ap is contained in the gauge-boson self- 
energies 

n z (o) n H ,(o) 

In the approximation considered, we write 

Ap = 3x t (l + 6 EW + 6 EW > QCD + 5 QCD ) ~ 3x t (l + p®x t )(l + h5^ D + h 2 6^' D ), (3) 
with 



x t = ITT" 1 ?: h = (4) 



\/2 2 a s 

^-m+, h — — . 

I671" Atx' 

a s being the QCD coupling constant. We have denoted by S EW the pure electroweak, 
by 8 EW >Q° D the mixed electroweak- QCD, and by 8® the pure QCD corrections. 
The two-loop electroweak correction p^ 2 \ due to virtual Higgs (ghost) effects, is small 
pW\ mH=0 = —0.74 for rriH ~ [0 but reaches a maximum as large as —11.57 at 
m H w 5.7m t @, |]. 

The one-loop correction to Ap was first calculated in [H. The two-loop QCD 
correction 5^ D has been calculated in ||. It proved to be rather large. If one takes 
m t as the top-quark pole mass, then 

= - fOr 2 + 3). (5) 



Therefore, it is essential to evaluate the next, three-loop correction, in view of the high 
precision of modern experiments. 



To evaluate Ap, the diagonal parts of the self-energies of the W and Z gauge bosons 



nn<?) = g^Mq 2 ) + q^tl a {q 2 ) (6) 

(a = W, Z) at q = are needed. Since at zero momentum and m ( ^ no infrared 
divergences appear in diagrams with only fermions and gluons, one may put q = from 
the very beginning. Contracting 11^(0) with g^, we obtain for IT Q (0) an expression 
containing only bubble integrals. At the one- and two-loop level these are quite simple 
and for arbitrary space-time dimension d = 4 — 2e can be written in terms of Euler's 
T function. Here we need only 

d% (m 2 )^ 2 _ rQg - d/2) 

d d k 1 d d k 2 ( m 2y+p +1 -d 



7r d {kl + m 2 ) a (k 2 + m 2 Y {{ki - k 2 ) 2 yi 
T{a + (3 + 1 -d) T(l - 7 ) r(/3 + 7 -f) T(a + j-j) 

r(a) r(p) r(a + p + 2 7 - d) 

d d h d d k 2 ( m 2^+/3+ 7 - a ! 

V d (k 2 ) a ((ki - k 2 ) 2 Y(k 2 + m 2 )"i 

T(a + (3 + 1 -d)T(a + (3- d/2) T(d/2 - a) T(d/2 - (3) 

~ r(a) r(/3) r( 7 ) r(d/2) ' 



(7) 



At the three-loop level 22 diagrams of the Z-boson self-energy and 29 diagrams of 
the H^-boson self-energy contribute to 5® CD . The integrals that appear here are much 
more complicated than at the one- and two- loop level. 

The rather complicated task of computing massive three-loop Feynman diagrams 
is accomplished by applying the method of recurrence relations |], |10[. This method 
allows us to relate various scalar Feynman integrals of the same prototype which differ 
by powers of their scalar propagators. As a result, by means of plain algebra, any 
diagram is reduced to a limited number of so-called master integrals. They need to be 
evaluated once and for all, and can then be used in any renormalizable quantum field 
theory. Some of the integrals that we need for the present three-loop calculation were 



considered in |10| . Here, however, more types of integrals are required. In addition 
to the master integrals evaluated in [10], two more nontrivial master integrals are 
encountered: 

d d h d d k 2 d d k 3 (m 2 ) 4 - 3d / 2 



[7r d / 2 r(l + e)f jfe?[(jfei - k 2 ) 2 + m 2 ] [(k 2 - k 3 ) 2 + m 2 ] [kl + m 2 ] 
1 15 65 135 81 „ s 

d d k X d d k 2 d d k 3 ( m 2)6-3d/2 
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[7r d / 2 r(l + e)} 3 kf(h - hykjKh - k 3 ) 2 + m 2 \\k 2 2 + m 2 \\{k 2 - k 3 ) 2 + m 
= 2C(3) - +D 3 + 0(e) (9) 

£ 



with 



So 



4 7T 

V3 Cl2( 3 ; 



0.260 434 137 632 162 098 955 729 . 



(10) 



We have not found a representation of D3 in terms of known transcendental numbers, 
though we do not exclude its existence. By means of the numerical method for the 
evaluation of Feynman diagrams proposed in [0 , D 3 can be calculated quite accurately. 
Here we give the first 22 digits, which is more than enough for a precise evaluation of 
P- 

D 3 = -3.027 009 493 987 652 019 786. (11) 



Calculations were mostly done using FORM 1.1 [|TTJ. All the diagrams were com- 
puted in the covariant gauge with an arbitrary gauge parameter. Performing charge 
and mass renormalization in the MS scheme, we got for the W-boson propagator the 
following expression: 



n£,'(o) = i2i,M^. 
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(12) 



Here n/ is the total number of quarks, / = ln(/i 2 /m 2 (//)), iritis) is the renormalized 
mass in the MS" scheme, and the constant 



Ba = 16 LL 



-ln 4 2--7r 2 ln 2 2- — 
3 3 180 



7T 



-1.762 800 087 073 770 086 (13) 



was defined in ||10|| . The result for the Z-boson propagator is 

3 5 
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(14) 



In the sum of the bare diagrams contributing to II a (0) (a = W,Z), as well as in the 
corresponding counterterms separately, the gauge parameter cancels, which is a partial 
check of our result. 

Substitution of (p~2|) and (111) into (0) gives the ultraviolet-finite expression 



S QCD _ r 2 
°(3),MS ~ ^ F 
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(15) 



If we perform mass renormalization in such a way that the renormalized mass is 
the pole mass m t , then we obtain 



cQCD 
0(3) 



c 2 
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(16) 



In the above formula I = ln(/i 2 /m 2 ) . Ultraviolet finiteness of (|15|) and (|Tj|) is an 
additional check of our result. Expression (1TB) can also be obtained from our result 



in the MS scheme by using the relation between rh t in the MS scheme and the pole 
mass mt 



13 



In the MS scheme for QCD [i.e. for the SU(3) gauge group with Ca = 3, Cp = 4/3] 
we get the following expression: 



xQcp 

MS 



1144 



16 

8 - — C(2) + 8/ ) h 



26459 16 „ 25064 
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(17) 



Substituting numerical values for all the constants and taking /z 2 = rhf with the 
minimally subtracted mass we obtain 

5®fL D = -0.061 511 928 430 2 a B -(0.221 937 307 314 4+0.030 043 860 323 8 n/) a 2 , (18) 

f = 6 turns into 



MS 

which at n 



xQCD 
MS 



-0.061 511 928 430 2 a s + 0.402 200 469 257 5 at. 



(19) 



The smallness of this correction in the MS scheme confirms expectations about higher 
order effects in electroweak parameters (see e.g. ||14|| ). 

With the definition of the renormalized mass m t as the pole mass of the top quark 
we obtain: 



fiQCD 
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Substituting numerical values for all the constants and putting fi 2 = mf we get 
5 qcd = _ .9io 338 291 586 9 a s -(2.564 571 412 664 2-0.180 981 195 7679 n f ) a 2 s , (21) 
and at nf = 6 we have 



13 



4 



11 



a 



~ -(20) 

7T 



fiQCD 



-0.910 338 291 586 9 a, - 1.478 684 237 779 7 at. 



[22) 



Both two- and three-loop QCD contributions are negative, and their effect is a 
screening of the bare mass splitting, so that only a reduced "effective" quantity enters 
the p parameter. 

Following the method of fastest apparent convergence 
loop correction into a rescaling of a s . With n./ = 6, tfju 
fj, ~ 0.2327 m t . When we apply the BLM procedure Wi 
in ( p0|) can be absorbed into the rescaling of a s 
same value was also obtained in [|17 



we can absorb our three- 
will be zero, if we take 
the n ^-dependent term 
if we choose p ~ 0.154 m t . The 
We conclude that the expression for the term 
proportional to n/, given in ]17|, agrees with the analytical result for this term given 
in (0). As an important result of our calculation, we stress the stability of S® CD : the 
usual perturbation theory, the FAC and BLM procedures give rather close results for 
5 QCD . Taking a s (m t ) = 0.1055, we have 5 QCD = -0.1125, -0.1159 and -0.1154 for 
the perturbation theory, FAC and BLM procedures, respectively. 

For our calculations we used the anticommuting 75. To check our result obtained 
with this prescription at least partially, we calculated S QCD again, using the regu- 
larization by dimensional reduction |TJj which keeps the algebra of 7 matrices four- 
dimensional. For propagator-type diagrams in the three-loop approximation the incon- 
sistency of this recipe is not yet revealed. As was expected, the result that we obtained 
in the regularization by dimensional reduction agrees with in the conventional 

dimensional regularization after the following recalculation of the coupling constant: 



h RD R = h(l+ -C A h 



(23) 



which just corresponds to a change in the renormalization scheme. The relation can 
be derived, for example, by equating invariant charges in these two schemes. 

Special attention was paid to the evaluation of the diagram with the axial anomaly. 
Only one such diagram contributes to Ap (namely, to 11^(0)) with two triangles and 
only the top quark running around. We evaluated it in the framework of the reg- 
ularization by dimensional reduction and using the prescription given in [[19]. The 



contribution from this diagram is finite, and the calculations yielded the same value in 
both approaches (although the evanescent parts were different). Our result also agrees 
with the one given in [p0|| . At n/ = 6 this contribution amounts to about 30% of the 
total three- loop correction (|22|). 

Some observables that are affected by our result are briefly mentioned here. One 
of them is the mass of the W boson as predicted from a, and Mz [21 



M 2 W 



pM. 



1 + 
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AAl 
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(24) 



where A Q 



ix a 



1/2 



37.2802(3) and Aa ~ 0.06 is the shift of the fine structure 



\V2G fl/ 

constant a due to photon vacuum polarization effects. The ellipsis stands for the 
non-leading remainder terms. Another physical quantity is the effective weak mixing 
parameter relevant to Z-resonance physics. It is given by 



sin 2 © 



M 2 W 
pM\ 



1 - 



\ PM, 



2 M-Aa 



+ 



(25) 



Numerical values illustrating how these observables are affected by various corrections 
are given in Table 1. QCD corrections were calculated using (p2]). 

Table 1: Percentage of various two- and three-loop QCD heavy-top corrections at 
a s (m t ) = 0.1055 and m t = 174 GeV. The value of a s was obtained by extrapolating 
a s (Mz=91.1895 GeV)=0.118 to the scale /i = m t with the aid of the three-loop (3 
function [2l| with n/ = 5. 



Observable 


Two-loop 
QCD 


Three-loop 
QCD 


Two-loop electroweak 
m H /m t = 1.5 5.7 10 


M w 
sin 2 9 


-0.065 
0.130 


-0.011 
0.022 


-0.002 -0.018 -0.025 -0.023 
0.003 0.036 0.051 0.046 



Table 1 demonstrates that the three-loop QCD correction is comparable with the 
two-loop electroweak correction for sizable Higgs masses (for mn/Tnt— 1.5 the former 
amounts to more than 60% of the latter). Thus, we conclude that it makes sense to 
evaluate subleading electroweak two-loop corrections to Ap (or another observable) 
only if the three-loop QCD corrections are taken into account as well. 
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